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Abstract 

The generation of large-scale magnetic fields is studied in dilaton electromagnetism in inflation- 
ary cosmology, taking into account the dilaton's evolution throughout inflation and reheating until 
it is stabilized with possible entropy production. It is shown that large-scale magnetic fields with 
observationally interesting strength at the present time could be generated if the conformal invari- 
ance of the Maxwell theory is broken through the coupling between the dilaton and electromagnetic 
fields in such a way that the resultant quantum fluctuations in the magnetic field have a nearly 
scale-invariant spectrum. If this condition is met, the amplitude of the generated magnetic field 
could be sufficiently large even in the case huge amount of entropy is produced with the dilution 
factor ~ 10 as the dilaton decays. 

PACS numbers: 98.80. Cq, 98.62.En 
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I. INTRODUCTION 



It is well known that magnetic fields are present on various scales in the Universe, from 
planets, stars, galaxies, to clusters of galaxies (for recent detailed reviews see [1-3]). The 
origin of the cosmic magnetic fields, however, is not well understood yet. Since they have 
direct influence not only on various astrophysical situations but also on the evolution of the 
Universe, their origin is one of the most important problems in modern cosmology. 

In galaxies of all types, magnetic fields with the field strength ~ 1CU 6 G, ordered on 
1 — lOkpc scale, have been detected 0, 0|. There is some evidence that they exist in galaxies 
at cosmological distances 0. Furthermore, in recent years magnetic fields in clusters of 
galaxies have been observed by means of the Faraday rotation measurements (RMs) of 
polarized electromagnetic radiation passing through an ionized medium 0. Unfortunately, 
however, RMs inform us of only the value of the product of the field strength along the 
line of sight and the coherence scale, and so we cannot know the strength without assuming 
the value of the coherence scale and vice versa. In general, the strength and the scale 
are estimated on 10~ 7 — 10~ 6 G and lOkpc— IMpc, respectively. It is very interesting and 
mysterious that magnetic fields in clusters of galaxies are as strong as galactic ones and that 
the coherence scale may be as large as ~Mpc. 

Some elaborated magnetohydrodynamical (MHD) mechanisms have been proposed to 
amplify very weak seed magnetic fields into the ~ 10 _6 G fields generally observed in galax- 
ies. These mechanisms, known as galactic dynamo 0], are based on the conversion of the 
kinetic energy of the turbulent motion of conductive interstellar medium into magnetic en- 
ergy. Galactic dynamo, however, is only an amplification mechanism, and so requires initial 
seed magnetic fields to feed on. Moreover, the effectiveness of the dynamo amplification 
mechanism in galaxies at high redshifts or clusters of galaxies is not well established. 

Scenarios for the origin of seed magnetic fields fall into two broad categories. One is 
astrophysical processes and the other is cosmological physical processes in the early Universe. 
The former, by and large, exploits the difference in mobility between electrons and ions. This 
difference can lead to electric currents and hence magnetic fields. The latter can also generate 
magnetic fields. Typically, magnetogenesis requires an out-of-thermal-equilibrium condition 
and a macroscopic parity violation. These conditions could have been naturally provided by 
the first-order cosmological electroweak phase transition (EWPT) [sj] or quark-hadron phase 
transition (QCDPT) [9( (see more references in the review j2j). The bubbles of new phase 
were formed in the old one and strong, though small-scale, turbulent motion is excited in 
the plasma. In standard model, however, EWPT is the second order, so that such bubbles 
cannot be formed. Furthermore, it has recently been shown by Durrer and Caprini ^(j that 
causally produced stochastic magnetic fields on large scales, e.g., during EWPT or even 
later, are much stronger suppressed than usually assumed. 

If the scale of cluster magnetic fields is as large as ~Mpc, it is likely that the origin of 
such a large-scale magnetic field is in physical processes in the early Universe rather than in 
astrophysical processes. From the two points, (1) There exists magnetic fields with the field 
strength ~ 10~ 6 G even in the objects where the effectiveness of the dynamo amplification 
mechanism is not well established, and (2) There is the possibility that the scale of cluster 
magnetic fields is as big as ~Mpc, it is conjectured that large-scale strong magnetic fields are 
produced in the early Universe and then are trapped in the plasma that collapsed to form 
galaxies and clusters of galaxies through adiabatic compression, or, in addition, secondary 
amplification mechanism such as galactic dynamo. 
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Since the conductivity of the Universe through most of its history is large, the magnetic 
field B evolves conserving magnetic flux as B oc a -2 , where a(t) is the scale factor. On the 
other hand, the average cosmic energy density p evolves as p oc a -3 in the matter dominated 
epoch. Hence B oc p 2 / 3 . The present ratio of interstellar medium density in galaxies p ga i to p 
and that of inter-cluster medium density in clusters of galaxies p cg are p ga i/p — 10 5 — 10 6 and 
p cg /p ~ 10 2 — 10 3 , respectively. Consequently, from these relations, we see that the required 
strength of the cosmic magnetic field at the structure formation, adiabatically rescaled to the 
present time, is 10~ 10 — 10 _9 G in order to explain the observed fields in galaxies B g8i \ ~ 1CT 6 G 
and clusters of galaxies B cg ~ 10 _7 G. On the other hand, in general, seed fields with the 
present strength 10~ 22 — 10~ 16 G is required for the galactic dynamo scenario. 

Although first-order cosmological phase transitions in the early Universe generate mag- 
netic fields, the comoving coherence length of the magnetic fields cannot be larger than the 
Hubble horizon at the phase transition, which is much smaller than Mpc today. Though the 
coherence length may grow due to MHD effects, this happens at the expense of the magnetic 
field strength. 

The most natural mechanism overcoming the large-coherence-scale problem is inflation in 
the early Universe (for a comprehensive introduction to inflation see Refs. 11, 3)- Turner 



and Widrow [13j (TW) first indicated that large-scale magnetic fields could be generated in 
the inflationary stage. Inflation naturally produces effects on very large scales, larger than 
Hubble horizon, starting from microphysical processes operating on a causally connected 
volume. If electromagnetic quantum fluctuations are amplified during inflation, they could 
appear today as large-scale static magnetic fields. This idea is based on the assumption that 
a given mode is excited quantum mechanically while it is subhorizon sized and then as it 
crosses outside the horizon "freezes in" as a classical fluctuation. However, there is a serious 
obstacle on the way of this nice scenario as argued below. 

It is well known that quantum fluctuations of massless scalar and tensor fields are very 
much amplified in the inflationary stage and create considerable density inhomogeneities 
13 | or relic gravitational waves [l5[. This is closely related to the fact that these fields 



are not conformally invariant even though they are massless. The amplification of the 
quantum fluctuations can be understood as particle production by an external gravitational 
field. Since the Friedmann-Robertson- Walker (FRW) metric usually considered is known 
to be conformally flat, the background gravitational field does not produce particles if the 
underlying theory is conformally invariant [16]. This is the case for photons since the classical 
electrodynamics is conformally invariant in the limit of vanishing masses of fermions. Hence 
electromagnetic waves could not be generated in cosmological background. If the origin of 
large-scale magnetic fields in clusters of galaxies is electromagnetic quantum fluctuations 
generated and amplified in the inflationary stage, the conformal invariance must have been 
broken at that time. Several breaking mechanisms have been proposed, which are mainly 
classified into the following three types. 

(1) A non-minimal coupling of electromagnetic fields to gravity: TW introduced the grav- 
itational couplings RA^A^, R pv A^A v , RF^F^/m 2 , etc, where R is the curvature scalar, 
A^ the electromagnetic potential, F^ the electromagnetic field-strength tensor, and m a 
constant with dimension of mass. The RA 2 terms could generate large-scale magnetic fields 
with interesting strength, but they also break gauge invariance by giving the photon an 
effective mass. In contrast, the RF 2 terms are theoretically more plausible, but the strength 
of the resultant magnetic fields is very weak. 

(2) A coupling of a scalar field to electromagnetic fields: TW first indicated the coupling 
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of an axion field, or that of a charged field which is not conformally coupled to gravity. After 
that many authors have studied more natural and effective couplings. 

Ratra 17( suggested the coupling of the inflation-driving scalar field (inflaton) (f) in the 
form e^F^F^ , and calculated the strength of large-scale magnetic fields in so-called a 
power-law inflation model induced by the exponential potential of the form e w< ^, where 
u and uj are constant parameters with dimension (mass) -1 . As a result, he found that 
present magnetic fields as large as 10~ 10 — 10~ 9 G could be generated. Recently Giovannini 
jl 8l | discussed the coupling of a massive scalar field ip m other than the inflaton in the form 
(ip m / Mp\)^ F^F^ , where £ is a constant parameter and Mp\ is the Planck mass. According to 
Giovannini, large-scale magnetic fields with the strength larger than the dynamo requirement 
could be generated. 

Garretson, Field, and Carroll analyzed the amplification of electromagnetic fluctuations 
by their coupling to a pseudo Goldstone boson (PGB) ip g in the form ip % F llv F tlu , where F^ v = 
1/2 e^ vpa F pa is the dual tensor of F^, and found that this coupling leads to exponential 
growth not for super-horizon modes but only for sub-horizon modes. Consequently, large- 
scale magnetic fields with interesting strength could not be generated |l9|. 

Magnetic fields due to a charged scalar field were considered in a special model in |2( 
(for more detailed review see |2l|). The authors found that stochastic currents could be 
generated during inflation due to production of charged scalar particles by the inflaton, and 
in turn, magnetic fields. Moreover, Davis et al. argued that the backreaction of the scalar 
field gives the gauge field an effective mass thus breaking the conformal invariance [22]. 
According to Davis et al., magnetic fields with the strength of order 10~ 24 G on a scale of 
lOOpc could be generated. 

(3) The conformal anomaly in the trace of the energy-momentum tensor induced by quan- 
tum corrections to Maxwell electrodynamics: It is known that the conformal anomaly, which 
is related to the triangle diagram connecting two photons to a graviton, breaks the confor- 
mal invariance by producing a nonvanishing trace of the energy-momentum tensor. Dolgov 
[23t | pointed out that such an effect may lead to strong electromagnetic fields amplification 
during inflation. According to Dolgov, however, magnetic fields with interesting strength 
might not be generated in realistic case, e.g., the model based on SU(5) gauge symmetry 
with three-fermion families. 

Incidentally, it has been indicated by Bertolami and Mota |24| that the conformal invari- 
ance might be broken actually due to the possibility of spontaneous breaking of the Lorentz 
invariance in the context of string theories, and that generated large-scale magnetic fields 
could be strong enough to explain the observed fields through adiabatic compression. 

In the light of the above various suggestions, it seems at present that the most natural 
and effective way of breaking the conformal invariance is to introduce the coupling of a scalar 
field to electromagnetic fields. In particular, Ratra's suggestion is attractive who claimed 
present magnetic fields as large as 10~ 10 — 10~ 9 G could be generated in his model, which 
would not require any dynamo amplification to account for the observed fields in galaxies 
and clusters of galaxies. 

In the present paper, in addition to the inflaton field we assume the existence of the 
dilaton field and introduce the coupling of it to electromagnetic fields. Such coupling is rea- 
sonable in the light of indications in higher- dimensional theories, e.g., string theories. Then 
we investigate the evolution of electromagnetic quantum fluctuations generated through 
the coupling, which breaks the conformal invariance of electrodynamics, and estimate the 
strength of large-scale magnetic fields at the present time. Particularly, we consider the 
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following two cases. One is the case the dilaton freezes at the end of inflation in the same 
way as Ratra [T3| just for comparison, and the other is the more realistic case that it still 
evolves after reheating and then decays into radiation with or without entropy production. 

Here we emphasize the following point. In Ratra's model, the inflaton and the dilaton are 
identified and power-law inflation is realized by introducing an exponential potential. There 
is no reason, however, why we should identify these fields. Furthermore, in the standard 
inflation models inflation is driven by the potential energy of the inflaton as it slowly rolls 
the potential hill. This slow roll over quasi-de Sitter stage is practically necessary to account 
for the nearly scale-invariant spectrum 1 of the primordial curvature perturbation out of the 
quantum fluctuations of the inflaton. 

The reset of this paper is organized as follows. In Sec. II we describe our model action 
and derive the equations of motion from it. In Sec. Ill we investigate the evolution of 
electromagnetic fields, and then estimate the strength of the large-scale magnetic fields at 
the present time in Sec. IV, where we assume the dilaton freezes at the end of inflation. On 
the other hand, in Sec. V, we consider the case the dilaton still evolves after reheating and 
then decays into radiation with or without entropy production. Although we consider the 
evolution of electromagnetic fields in slow-roll exponential inflation models in Sees. II— V, 
for comparison we discuss it in power-law inflation models in Sec. VI keeping the recent 
WMAP data in mind. Finally, Sec. VII is devoted to discussion and conclusion. 

We use units in which &b = c = H = 1 and denote the gravitational constant 8kG by k 2 
so that k 2 = 8-7r/M P1 2 where M P1 = G~ x l 2 = 1.2 x 10 19 GeV is the Planck mass. Moreover, in 
terms of electromagnetism we adopt Heaviside-Lorentz units. The suffixes 'i', '1', 'R', and 
'0' represent the quantities at the initial time t> u the time when a given mode first crosses the 
horizon during inflation t\, the end of inflation (namely, the instantaneous reheating stage) 
tR, and the present time to, respectively. 



II. MODEL 
A. Action 

We introduce two scalar fields, the inflaton field 0, and the dilaton $. Moreover, we 
introduce the coupling of the dilaton to electromagnetic fields. Our model action is given as 
follows. 

S = J d 4 X\^g [ £inflaton + ^dilaton + ^EM ] , (1) 

Anflaton = " ^T<^<R0 - U[<j>\ , (2) 

^dilaton = -~<r0 M $0„$-V[$], (3) 

Cem = -\f{®)F^, (4) 



1 The spectral index is estimated as 0.99 ± 0.04 by using the first year Wilkinson Microwave Anisotropy 
Probe (WMAP) data only [2^, where the errors are the 68% confidence interval. 
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/($) = exp(A^), (5) 
V[<£>] = Vexp(-\K$), (6) 

where g is the determinant of the metric tensor g^, U[<p] and V[$] are the inflaton and 
dilaton potentials, V is a constant, and / is the coupling between the dilaton and electro- 
magnetic fields with A and A ( > ) being dimensionless constants. The form of the coupling 
between the dilaton and electromagnetic fields in Eq. and that of the dilaton potential 
in Eq. © can be motivated by higher-dimensional theories reduced to four dimensions. 

We assume the spatially flat Friedmann-Robertson- Walker (FRW) space-time with the 
metric 

ds 2 = g (W dx tl dx v = -dt 2 + a 2 (t)dx 2 , (7) 

where a(t) is the scale factor. In terms of the U(l) gauge field A^, the electromagnetic 
field-strength tensor is given by 

= d„A v - d u A^. (8) 
Before going on, we state the framework we adopt in this paper. 

(1) During slow- roll inflation the cosmic energy density is dominated by U[<f>] and the 
energy density of the dilaton is negligible. 

(2) The Universe is reheated immediately after inflation at t — £r. See e.g. |26| for an 
efficient mechanism of reheating. 

(3) The conductivity of the Universe a c is negligibly small during inflation, because there 
are few charged particles at that time. After reheating a number of charged particles are 
produced, so that the conductivity immediately jumps to a large value: cr c ^> H (t > t-&), 
here H is the Hubble parameter. This assumption is justified by a microphysical analysis 

0. 

(4) We consider both the case the dilaton is frozen at t — t& as in Ratra's model [Tj| and 
that it continues to evolve after reheating until it decays with or without entropy production 
and gets stabilized at a potential minimum. 

(5) The value of the coupling / between the dilaton and electromagnetic fields is set to 
unity when the dilaton gets stabilized so that the standard Maxwell theory is recovered: 

/ = !■ 



B. Equations of motion 



From the above action in Eq. the equations of motion for the inflaton, the dilaton, 
and electromagnetic fields can be derived as follows. 



1 -a w=5riu) + 



-g 



1 -A W=5,TM) ■ dV[i] 



-g 



,d, (V=?/($)F^) = 0. 



0. 



4 d$ ^ ' 



(9) 
(10) 

(11) 
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Since we are interested in the specific case where the background space-time is inflating, 
we assume that the spatial derivatives of <fi and $ are negligible compared to the other 
terms (if this is not the case at the beginning of inflation, any spatial inhomogeneities will 
quickly be inflated away and this assumption will quickly become very accurate). Hence the 
equations of motion for the background homogeneous scalar fields read 

+ 3if0+^M = O, (12) 



$ + 3#$ + = 0, (13) 



together with the background Friedmann equation 



2 2 

a \ k 



= [a) = T( p + + P * ) ' (M) 
P* = \ft + U[<f>], (15) 

p* = ^> 2 + m (i6) 

where a dot denotes a time derivative. Here p^ and p$ are the energy density of the inflaton 
and that of the dilaton. Since we have p^ ^> p$ by assumption, during inflation H reads 

H 2 « = i/inf 2 , (17) 

where H in f is the Hubble constant in the inflationary stage. 

We consider the evolution of the gauge field in this background. Its equation of motion 
in the Coulomb gauge, A (t, x) = and djA^it, x) = 0, becomes 



Mt, x)+\H + J -)A i {t,x)- -d^Mt, x) = 0. (18) 



III. EVOLUTION OF MAGNETIC FIELDS 

In this section, we investigate the evolution of the vector potential and then consider that 
of the electric and magnetic fields. 



A. Evolution of vector potential 

To begin with, we shall quantize the vector potential Ai(t,x). It follows from the elec- 
tromagnetic part of our model Lagrangian in Eq. (@J) that the canonical momenta conjugate 
to the electromagnetic potential A M (i, x) are given by 

vr = 0, 7r i = /($)a(t)A i (t,aj). (19) 
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We impose the canonical commutation relation between Ai(t,x) and 7Tj(t,x): 



[Ai(t,x),7ij(t,y)] = i 



d k r ih-(m-y) fx _ tih. 
(27T) 3 V 3 V 



(20) 



where k is co moving wave number, and k denotes its amplitude |fe|. From this relation we 
obtain the expression for Ai(t, x) as 



Ai(t,x) 



d 3 k 



b(k)Ai(t, k)e ik ' x + P{k)A*{t, k) e - ik 



(2tt) 3/2 

where b(k) and b\k) are the annihilation and creation operators which satisfy 



(21) 



b(k),b\k') =5 3 {k~k'), b(k),b(k') = P(k),P(k') 



0. 



(22) 



From now on we choose the x 1 axis to lie along the spatial momentum direction k and 
denote the transverse directions x 1 with I = 2, 3. From Eq. ([18)1 we find that the Fourier 
modes of the vector potential Ai(t, k) satisfy the following equation: 



./ 



Ar(t, k) + \ H M + j ) Aj(t, k) + — k) = 0, 



(23) 



and that the normalization condition for Ai(t, k) reads 



Ai(t, k)A*(t, k) - Aj(t, k)Ai*(t, k) = — 8,, 

fa 



(24) 



For convenience in finding the solutions of Eq. (|23jl . we introduce the following approxi- 
mate form as the expression of /. 



/($) = /[$(t)] = /[$(o(t))]= fa 



13-1 



(25) 



where / is a constant and (3 a parameter whose time-dependence is weak as will be seen in 
the next subsection. Using Eq. (J23j) . we find 



/ 

#inf + J = /3-Hinf ■ 



(26) 



It follows from Eq. (j2Tj) that Eq. (J2*3*j) is rewritten to the following form by replacing the 
independent variable t to rj. 



(fA^k.r]) (1- (3\ dAi(k, 77) 



drf 



dr) 



+ k 2 A I (k,i 1 ) = 0, 



(27) 



where 77 = J dt/a(t) is conformal time. During inflation 77 = —l/(aH m {). If we regard /3 as 
a constant, the solution is given by 

Aj{k,r]) = Cj+ikX-H^Y^Hfyi-kri) + Cj^-H^f^H^-krj), (28) 
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where is an u-th order Hankel function of type n (n = 1,2), and Ci + (k) and Ci-(k) 
are constants which satisfy 

\C I+ {k)\ 2 -\CUk)\ 2 = (29) 

4-Hinf J 

We shall choose 



Ci+{k) = e^'\ CUk) = 0, (30) 

so that the vacuum reduces to the one in Minkowski space-time at the short-wavelength 
limit: 

kr] ^oo; A I (k,7 1 ) -> -^e^. (31) 



We therefore obtain 

Being interested in large-scale magnetic fields, we investigate the behavior of this solution 
in the large- wavelength limit. Expanding the Hankel function in Eq. (}3*2"j) and taking the 
first leading order in fc/(aifi n f), we obtain 



-P/2 



e i(f3-l)n/4 



AAk,a) = 7P'\ ( £) a' 1 ' 2 ( — Y 

A ' ' ]j4nH ini f(a) \2 ) \aH ini J 

oc k~ p/2 a , for (3 > 0, (33) 



and 



P/2 



e i(3-/3>/4 



4irH in{ f(a) \ 2 J \aH inf/ 

oc A^ /2 a"' 3 , for /? < 0. (34) 

The large-scale expansion bifurcates at /3 = 0. This is the reason for the two different 
expressions, Eqs. (jHSJ) and (pMJ) . From Eqs. (|3l5|) and (|HIJ), we see that the large-scale vector 
potential is time-independent for {3 > and evolves like a _/3 for /3 < 0. Furthermore the mean 
square of the large-scale vector potential in the position space is ~ k 3 \Aj(k, a)\ 2 oc AH^ -1-3 
on a comoving scale r = 2n/k, so the root-mean-square (rms) has a scale- invariant spectrum 
when 1/31 = 3. 



B. Parameter /3 

In the previous subsection, for convenience we have introduced the form (J23j) as the 
expression of /. In practice, however, it is a function of $ as seen in (jSJ). We therefore 
investigate the expression of parameter (3 by the comparison of Eq. © with Eq. (|23j) . 
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To begin with, we consider the evolution of the scale factor a(t) and the dilaton $(£) in 
the inflationary stage. The scale factor a(t) is given by 



a(t) = aiexp[H inf (t - h)}, 



(35) 



where a± is the scale factor at the time t\ when a given comoving wavelength 27r//c of the 
vector potential first crosses outside the horizon during inflation, k/(aiH- m {) = 1. In order 
to obtain the analytic solution of Eq. (JT3j) we apply slow-roll approximation to the dilaton, 
that is, 



<1, 



(36) 



and then Eq. (fT^J) is reduced to 



3#inf$ + 



0. 



(37) 



The solution of this equation is given by 



$ = — In 

Xk 



(Xk) 2 V 



3H; 



inf 



(t-t R ) +exp(A/€$ R ) 



(38) 



where $r ( < ) is the value of $ at the end of inflation. In the case the dilaton is frozen at 
t — tji, we choose $r = so that / = 1 at that time. 

Next, we investigate the evolution of ///. Using Eqs. (j5J), (fTTj) and (|38|). we find 



/ 



\\k 2 V[$] 



3H; 



XXH inf w, 



inf 



where u> is defined as 



w = 



P4> ~ ' 



Since we have ^> p$ by assumption, w < 1. Comparing Eq. (j2T?j) with Eq. 

P « 1 + AAw = 1 + Xe, 

where X and e are defined as 

e = A 2 u>, 



(39) 

(40) 

, we find 
(41) 

(42) 
(43) 



respectively. 
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Though U[(p} is approximately constant in slow-roll exponential inflation, V[<&] changes 
gradually. As a result it follows from Eqs. ©, P5|). P0jl. and (j4Hj) that the ratio of iu at 
t = £r to that at t = ti is 

™(*r) . , ^[$(*r)] 



«>(*i) 



, (Ak) 2 ^[$r] 



3i?inf' 



#inf(^R — ^i) 



(44) 



where c(£r) = A 2 u>(£r). Since -ffi n f(tR — 1{) = N(ti — >■ £r) is about 50, where N(U — > £r) is 
the number of e-folds during the period from i = t\ to i = £r, the ratio of w(£r) to is 
about half in the case e(tn) ~ 1/100. If e 1, the variation in w in the inflationary stage 
is small, and then we can regard (3 as approximately constant. 

The slow- roll condition to the dilaton, Eq. (|36|). is equivalent to the following relation. 



< 1 



e< 1. 



(45) 



If we assume A ~ 0(1), the relation e 1 is satisfied during inflation because of w <C 1. 



C. Evolution of electric and magnetic fields 

We consider the evolution of electric and magnetic fields. The proper electric and mag- 
netic fields are given by 

Ei pmpei {t, x) = a- l Ei{t, x) = -a^Aiit, x), (46) 
Br peT {t, x) = aT l B^t, x) = a- 2 e ijk djA k (t, x), (47) 

where Ei(t,x) and Bi(t,x) are the comoving electric and magnetic fields, and is the 
totally antisymmetric tensor (e 12 3 = 1). 

From Eqs. (JH2J), (j4fijl . and (|4Tj) we find the Fourier components of the comoving electric 
and magnetic fields in the inflationary stage: 

s ' (W ^0"-" , *.(i)' ww ' <«> 

™ - o e,( " +iw4 <49) 

We consider the case the dilaton freezes at the end of inflation and after instantaneous 
reheating the conductivity of the Universe a c jumps to a value much larger than the Hubble 
parameter at that time. The evolutionary equation of the vector potential for an electrically 
conducting plasma is given by 

Ai{t, x)+(^ + a c ^j Ai(t, x) - -^OjdjAit, x) = 0. (50) 
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The joining conditions at the transition from the inflationary stage (INF) to the radiation- 
dominated one (RD) at t = t R are [17| 

E+ RD \t R ,x) = exp(-a c t R )E+ mF \t R ,x), (51) 

B^(t R ,x) = B^ F \t R ,x). (52) 

From these joining conditions, for a large enough conductivity at the instantaneous reheating 
stage, we see that the electric fields accelerate charged particles and dissipate. In fact, we 
solve Eq. (jBTIj) in the large-conductivity limit: a c 3> H, so that electric fields vanish and the 
proper magnetic fields evolve in proportion to a~ 2 (t) in the radiation-dominated stage and 
the subsequent matter-dominated stage (t > t R ) |l7| |. 

It follows from Eq. (}4T?j) that the Fourier components of the proper magnetic fields in the 
inflationary stage are given by 

Improper ^|2 _ Q -2 Q ) |2 

= (i^d/w) (a) Q) (^) H W (^fc) ' (53) 

Since we are interested in the scales much larger than the Hubble radius in the inflationary 
stage, we expand Eq. (J53j) in the large-scale limit. Taking account of the above fact that 
the proper magnetic fields evolve in proportion to a~ 2 (t) (t > t R ), we find that the Fourier 
components of the large-scale proper magnetic fields are expressed as 

| 5/ p-p- ( ^)| 2 = ^r 2 (tj f -\a R ) 



OR-f^inf/ V a R-^inf/ \ a 



kYfV 2 



for (3 > 0, (54) 



and 



9-(/ 3 + 2 ) / R\ 

\ Bl ^(t,k)\ 2 = —^r 2 f-|J r\a R ) 



for /3 < 0, (55) 



where we have only recorded the first leading term. 

Finally, the energy density of the large-scale magnetic fields in Fourier space is given by 

p B (t,k) = \B I ^(t,k)\ 2 f(a). (56) 

Multiplying ps(t,k) by phase-space density: 47ik 3 /(2ii) 3 , we obtain the energy density of 
the large-scale magnetic fields in the position space 

k 3 

p B (L,t) = —\B^(t,k)\ 2 f(a), (57) 
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on a comoving scale L = 2ir/k. Note that the expressions in Eqs. (|54|) and (jo3j) are for 
one transverse component of the magnetic fields and the total magnetic field contribution 
is twice larger. Using Eqs. (p>4"j). (fo3j). and (joTj) . we find that the energy density of the 
large-scale magnetic fields at the present time t is given by 



PB{L,t 



21/31-3 



TP' 



r 2 ( ^ ) HJ" f ^ 

a 



k 



inf 



-1/31+5 



(58) 



Hence the large-scale magnetic fields have a scale-invariant spectrum when 
now on we shall take the present scale factor ao = 1. 



5. From 



D. Consistency 

During inflation the energy density of electric and magnetic fields should be smaller than 
that of the dilaton so that the evolution of the dilaton is governed by its (classical) potential. 

We define the ratio of the energy density of the electric and magnetic fields to that of the 
dilaton as follows. 



T(L,t) = — [p B (L,t)+p E (L,t)} 
P<s> 



p E (L,t) 



k^ 

2^ 



£/ proper (U)| 2 /(a), 



(59) 
(60) 



where Pe{L, t) is the energy density of the electric fields on comoving scale L = 2ir /k. During 
inflation T must be smaller than 1, which we call the consistency condition. 
Using Eqs. (H7J), (gHJ), flU, flUJ), and (jUD we find 



T ss — 



inf 



3w \M Pl 



k 



inf 



H 



(i) 



k 



inf 



H 



(2) 



PI 2 \aH; 



inf 



(i) 



^Z 2 - 1 V aft 



k 



'inf 



H 



(2) 

v 2 - 1 V am 



k 



'inf 



(61) 



where the approximate equality follows from the ratio of the energy density of the dilaton 
to that of the inflaton in the inflationary stage: p^/ptf, ~ w, see Eqs. ffTtjj) and (pfUj) . 

Here we note that the upper limit on H m { is determined by the observation of the 
anisotropy of cosmic microwave background (CMB) radiation. Using the WMAP data on 
temperature fluctuation (Hf, which is consistent with the Cosmic Background Explorer 



(COBE) data, we can obtain a constraint on if in f from tensor perturbation [15, 27 



-f^inf 



< 2 x 1(T 5 . 



(62) 



pi 



From this relation we find that the upper limit on H m f is 2.4 x 10 14 GeV. 

When we consider a given scale 2n/k, the value of k/(aH in {) decreases as the Universe 
evolves. Evaluating Eq. (JBTj) at the horizon crossing k/(aH\ n {) = 1, we find T < 10 



-10 



w 
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and so the consistency condition is satisfied. In the long wavelength regime, expanding Eq. 
(|61|) and taking the first leading order in k/(aH- m t), we find 



T 



1 



3ir 2 w \Mpi 



X 



(3 



T 



+ 2 ^- 2 r 2 1 - - i 

2 1 I 2 



1 ( Hi ni \ 2 ( aH mi ^ 



aH-, 



inf 



-2 



for > 2, 



(63) 



3vr 2 w V Mi 



'pi 



x 



and 



T 



2 ^r 2 ( ^ + 2-^- 2 )r 2 + i 



1 {H mi V(aH^-^ 



inf 



-2(/?-l) 



for < /3 < 2, (64) 



37T 2 U> VMpi 







2 -/3r 2 f-|j +2-^- 2 )r 2 ( + i 



inf 



for /3 < 0, (65) 



where the first term in the square parentheses is the magnetic contribution and the second 
term the electric one. T is smaller on larger scales for —3 < (3 < 5. 

IV. ESTIMATION OF PRESENT LARGE-SCALE MAGNETIC FIELDS 

In this section, we estimate the present strength of the large-scale magnetic fields. Since 
we consider the case the dilaton is frozen at t = £r, we choose $r = so that / = 1 at that 
time. 

A. estimation of present large-scale magnetic fields 

To estimate the present large-scale magnetic fields from Eq. ()58|) we use the following 
values 13 ■ 



H = 70h 70 kms" 1 Mpc" 1 w 2.26/i 70 x HT 18 s -1 



7T 



30 

iV = 45 + In 



g*T R 4 (0„«2OO), 



or 



[Mpc 

g* \V3 t r 



[Mpc]; 1 



3.91 



[30/(vr 2 ^)] 1/12 p// 4 
| 10 38 / 3 [GeV] 

3.7T R 



'--yd 



2.35 x 10- 13 [GeV] 



(T 70 ^2.73K), 



(66) 
(67) 

(68) 
(69) 
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where H is the present Hubble parameter (throughout this paper we use h 70 = 1.00 |2f 
g* is the total number of degrees of freedom for relativistic particles at the reheating epoch, 
Tr is reheating temperature, N is the number of e-folds between the time t\ and the end of 
inflation £r, and T 70 is the present temperature of CMB radiation. 

Applying Eqs. (fo% Jl -(ffi3 |) and k/(a R H inf ) = exp(-iV) to Eq. we find 



ML, t ) K \ B ( L >to)\ 2 « (VtfinfMp^'" 5 | T 70 y^ ) . (70) 



From this relation we see that stronger magnetic fields could be generated in the case of 
a large H ini and a large \/3\. The reason is as follows. From Eqs. (J5Hj) and ()43|). we find 
$(t ; ) = (A/t) _1 ln[l - eifinf(*R - ti)]. Furthermore, from /// = ((3 — l)H mi , we find that 
the rate of the change of / is larger in the case of a large \/3\ and a large H in f. Hence the 
conformal invariance of the Maxwell theory is broken to a larger extent in the case of a large 
.f/inf and a large \/3\. Practically, for (3 > 0, the initial amplitude of electromagnetic quantum 
fluctuation becomes larger. On the other hand, for (3 < 0, although the initial amplitude 
smaller, the rate of the amplification is very high. We will consider the detailed values of 
the present magnetic fields in Sec. IV C. 



B. Upper limits on cosmological magnetic fields 

Upper limits on cosmological mag netic fields come from the following three sources (see 
more detailed explanations in |3l l29l|). 

(1) CMB anisotropy measurements: Homogeneous magnetic fields during the time of 
decoupling whose scales are larger than the horizon at that time cause the Universe to expand 
at different rates in different directions. Since anisotropic expansion of this type distorts 
CMB, measurements of CMB angular power spectrum impose limits on the cosmological 
magnetic fields. Barrow, Ferreira, and Silk [3fJ carried out a statistical analysis based 
on the 4-year COBE data for angular anisotropy and derived the following limit on the 
primordial magnetic fields that are coherent on scale larger than the present horizon. 

a£Lc < 5 x 10- 9 G. (71) 

Incidentally, Caprini, Durrer, and Kahniashvili j3lj have recently investigated the effect of 
gravity waves induced by a possible helicity-component of a primordial magnetic field on 
CMB temperature anisotropies and polarization. According to them, the effect could be 
sufficiently large to be observable if the spectrum of the primordial magnetic field is close 
to scale invariant and if its helical component is stronger than ~ 10~ 10 G. Moreover, it has 
also been argued that if tangled magnetic fields of > 3 x 10~ 9 G exist on cosmological scales, 
their imprint on CMB anisotropy and polarization maybe detectable [32l |. 

(2) Big Bang Nucleosynthesis (BBN): Magnetic fields that existed during the BBN epoch 
would affect the expansion rate, reaction rates, and electron density. Taking all these effects 
into account in calculation of the element abundances, and then comparing the results with 
the observed abundances, one can set limits on the strength of the magnetic fields. The 
limits on homogeneous magnetic fields on the BBN horizon size ~ 1.4 x 10- 4 /i 70 -1 Mpc are 
less than 10~ 6 G in terms of today's values |33^ . 

(3) Rotation Measure (RM) observations: RM data for high-redshift sources can be used 
to constrain the large-scale magnetic fields. For example, Vallee [34| tested for an RM 
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dipole in a sample of 309 galaxies and quasars. The galaxies in this sample extended to 
z ~ 3.6 though most of the objects were at z < 2. Vallee derived an upper limit of 
6 x 10~ 10 (n e0 /10~'cm~ 3 ) 1 G, where n e0 is the present mean density of thermal electrons, on 
the strength of uniform component of a cosmological magnetic field. Note that the average 
baryon density is estimated as = (2.7 ± 0.1) x 10~ 7 cm -3 [25[ for comparison. 



C. Results 

We show the results of the present large-scale magnetic fields calculated from Eq. (|58|). As 
described in Sec. I, primordial magnetic fields with the present strength 10 -10 — 10~ 9 G are 
required to explain the observed fields in galaxies and clusters of galaxies through adiabatic 
compression. On the other hand, for galactic dynamo scenario, the fields with the present 
strength 10~ 22 — 10~ 16 G are required. 

Here we define Q(L, t R ) as the ratio of the energy density of the large-scale electric fields 
to that of the magnetic counterpart at the end of inflation tn, 

Q(Lt )= |£/ proper (-MR)| 2 (72) 
' \B! pmpcT {L,t R )\ 2 ' V } 



From Eqs. flU - fU) , flUJ) and (J?2J), we find 

r 2 (/3/2-i) r _ 2/ 



©(4 *r) = jZhjq J K- 2 (L), for > 2, (73) 



and 



where 



Q{LM = V \2% 2 ,2) 1) k ~ 2{P ~ 1){L) ' for < /? < 2, (74) 



Q(L,t R ) = T \l 3/ l + 1) K 2 (L), for/5<0, (75) 



V\-(5/2) 



K<L) = &[mZg^\(L\. (76) 



\ 10 38 / 3 [GeV] JVtMpc 



While the energy density of magnetic fields is dominant for (3 > 1, that of the electric fields 
is dominant for (3 < 1 on large scales. Particularly, in the case (5 < 0, the energy density of 
the large-scale electric fields is much larger than that of the magnetic counterpart. When 
strong enough magnetic field is generated, the accompanying electric field is so strong that 
the consistency condition is not satisfied. The reason is as follows. As noted in Sec. Ill A, 
for (3 > 0, the vector potential is time-independent in the large-scale limit in the inflationary 
stage, conversely, for /3 < 0, it evolves like a _/3 . Hence the electric fields are not generated 
in the former case, but generated in the latter case. Consequently, stronger magnetic fields 
could be generated in the case (3 > 1. 

Figures 1 and 2 depict the curves in the H in f — (3 parameter space on which the present 
magnetic fields on IMpc scale with each strength could be generated. The former is for 
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(3 > and the latter for (3 < 0. In Fig. 1, the shaded area illustrates the excluded region 
from the observation of the anisotropy of CMB. This region is decided by Eq. ()7ip . From 
Eqs. (J63J) and we find that the consistency condition is satisfied in all the area. On the 
other hand, in Fig. 2, the shaded area corresponds to T > 1 and illustrates the excluded 
region from the consistency condition decided by Eq. (|65|) . The excluded region from CMB 
is included in this region. From these figures we can find the following results. For (3 > 0, 
the magnetic fields on IMpc scale with the strength larger than 10 _10 G could be generated 
in the case H m { > 1.2 x 10 6 GeV. Furthermore, the magnetic fields strong enough for the 
galactic dynamo scenario (10~ 22 — 10~ 16 G) could be generated in the wide region of the 
parameter space. On the other hand, for (3 < 0, the maximum strength of the fields are 
5 x 10 _19 G at most. Parenthetically we note that in all figures we depict 10~ 22 G contours 
with solid lines to emphasize the critical value for the galactic dynamo scenario. As we can 
see from Eq. (J7UJ), when (3 > 5, the larger-scale field is the stronger, so that its strength is 
constrained from the observation of the anisotropy of CMB, on the other hand, when (3 < 5, 
the stronger is the smaller-scale field, so that its strength is constrained from the predictions 
of light element abundances from BBN. The former is more stringent than the latter, and 
so all results shown in Figs. 1 and 2 satisfy the limit imposed by the latter. Here we show a 
few characteristic examples. When (3 = 1, i.e., in the ordinary Maxwell theory, the present 
strength of the magnetic field on IMpc scale is 2.1 x 10~ 58 G. This value is independent 
of -f/inf- Moreover, when (3 = 5 (the magnetic field has a scale-invariant spectrum), the 
maximum value of the field is 1.8 x 10 _11 G in the case H m ^ = 2.4 x 10 14 GeV. 

Finally, we note the following point. As discussed in Sec. Ill B, we have assumed e <C 1 
and regarded (3 as approximately constant, and analytically investigated the evolution of 
the electric and magnetic fields. This approximate analysis is proper. In fact we have 
taken account of the time-dependence of (3 in the case w{t^) <C 1 and A ~ 0(1)\ P ~ 
l + Xw{t R )/ (l + w{t R )H m{ (t-t R )), which is derived by using Eqs. ©, (001), and flUJ), 
and numerically solved Eq. (|2T|) in the inflationary stage. As a result we have confirmed 
that the numerical results almost agrees with the analytic ones. 

V. DILATON DECAY AFTER REHEATING 

So far we considered the case the dilaton field freezes at the end of inflation and so 
the value of the coupling / between the dilaton and electromagnetic fields is set to unity. 
In practice, however, it is expected that the dilaton continues its evolution along with 
the exponential potential even after reheating but is finally stabilized when it feels other 
contributions to its potential from, say, gaugino condensation that generates a potential 
minimum. As it reaches there, the dilaton starts oscillation with mass m and finally decays 
into radiation. Then we consider its potential minimum is located at $ = and / = 1 there, 
and so the field amplitude evolves from $r( < ) to nearby zero along with the exponential 
potential. 

During the coherent oscillation the energy density of the dilaton p$ evolves as a~ 3 (t), 
so that it decreases more slowly than that of the radiation produced by the inflaton p r ^ nf \ 
If the Universe is radiation-dominated until the dilaton decays, the entropy per comoving 
volume remains practically constant. If p$ becomes dominant over p^ mi \ significant amount 
of entropy is produced, which dilutes the energy density of magnetic fields. 
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A. Dilaton decay 



We regard the time t osc ~ m _1 as the epoch when the coherent oscillations commence. 
When t > £r, the Universe is radiation-dominated, and so we shall set a[t) = a^{t/t-p) 1 / 2 . 

Before the field oscillation regime (£r < t < t osc ), the dilaton evolves with the exponential 
potential (6). It is known that, if the dilaton evolves with the exponential potential for a 
sufficiently long time, the dilaton enters a scaling regime as it evolves down its potential (35^ . 
In this scaling regime the friction term from the expansion of the Universe in the equation 
of motion balances the potential force allowing it to enter this scaling era. If the Universe 
is radiation-dominated in this regime, therefore, the evolution of the dilaton potential is the 
same as that of the background radiation produced by the inflaton, which is proportional 

to r 2 . 

Then we have numerically solved the equation of motion for the dilaton along with the 
exponential potential, and found that the dilaton evolves near to zero before it enters a 
scaling era because the field amplitude at the end of inflation |<3>r| is constrained to be 
relatively small. The reason is as follows. When £r < t < t osc , the energy density of 
magnetic fields is enhanced through the coupling with the evolving dilaton. The energy 
density of magnetic fields on all scales should remain much smaller than that of the dilaton 
so that (4) should not affect the evolution of the dilaton. As we will show in Eq. (19 1J) in 
the next subsection, the energy density of magnetic fields is more enhanced in the case of a 
large value of Ak|$rJ. Hence the upper limit on Ak|$r| is determined from this condition. 
We have therefore numerically solved the equation of motion for the dilaton along with 
the exponential potential and investigated the region of Ak|$rJ in which this condition is 
satisfied. As a result we have found that the value of |$r| must be relatively small as 
I^rI ~ l//c in the case A ~ From now on we will deal with Ak|$rJ as a model 

parameter, and show the detailed region of this parameter in Figs. 3 and 4 in Sec. V C. 

Consequently, even if the energy density of radiation pf°^ is dominant over that of the 
dilaton p$ at the instantaneous reheating stage, that is, V"exp(— Ak$r,)/p0 «w<1, p? ni ^ 
becomes comparable to p$ around the field oscillation regime (t ~ t OBC ). Since p$ does not 
become dominant over pr inf ^ in £r < t < t osc , however, it is the field oscillation regime that 
significant amount of entropy could be produced. 

We consider the evolution of p$ and the energy density of the radiation produced by the 
dilaton decay p r ( dil ) in the epoch of the coherent oscillations (t > t osc ). The equations for 
them are given as follows 2 • 



P* 



3- + T$ ) p$, 

a 



(77) 



P 



(dii) 



4-p( dil ) + r 



$p$, 



(78) 



where r$ is the decay width of $. The solutions of these equations are given by 

l -3 

exp [—T 9 (t - tosc)] , 



p$ — p<s>{t osc ) 

c 



a(t) 



a(t a 



(dii) 



a(t) 
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_o(tosc) _ 




_Ct(t osc ) _ 



exp [-r$(r - W)] dr. 



(79) 
(80) 
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Here we estimate p<$>{t = t OBC ) ~ V because V is the value of the exponential potential in 
the form © at $ = 0, and we expect the contribution from stabilization mechanism has the 
same order of magnitude. 



B. Entropy production 



If the Universe is radiation-dominated until the dilaton decays t$ ~ T$ , we find the 
following relation. 



or 



a(t Qf 



f > rv 1/2 *R- 2 m- 3/2 

P<t> 



V 



> 



M 



pi 



m 



inf 



m 



(81) 



where the last relation is obtained by using T$ ~ m(m/Mp\) 2 and £r « 1 / (2iJ iri f ) . In this 
case the entropy per comoving volume remains constant. Hence the necessary condition of 
entropy production is 



M 



pi 



m 



2H-, 



inf 



m 



(82) 



Here we consider the case p r ( inf ) becomes equal to p$ at t = t c in the epoch of the coherent 
oscillations (t osc < £ c < £$): 



Pr (ini \t c )=p^(t c ) 

From this equation we find 



a(t c ) 



P$(t c 



a(t c 
_a(t os 



(83) 



' < ~ ( J £r t osc 



'inf 



(84) 



where the first approximate equality follows from p<p(t osc ) ~ V, and the second one from 
£ osc ~ m^ 1 and £r « l/(2f/~i n f). When £ > £ c , p$ is dominant over p r ^ mi \ so that the Universe 
is matter-dominated. Thus we shall set a(t) = a c (£/£ c ) 2 / 3 = a^t^i' 1 ^ 2 t c ~ 1 ^ 6 t 2 ^, where the 
second equality follows from the joining condition of a(t) at £ = t R . 

We investigate the entropy per comoving volume after the dilaton decay. In general, the 
entropy per comoving volume is represented as S = a 3 (p + p)/T, where p, p and T are the 
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equilibrium energy density, pressure, and temperature, respectively. It is given by 



S 4/3 = Sc 4/3 + ^*(*cK 4 



45 



n l/3 



r„ 



t r 



a(r) 



exp [— r$(r — t c ) ] dr (85) 



/■oo 

1 + rV c ~ 2/3 / (u + t c ) 2/3 exp (-r*u) dw 



i + r ( - ) (t c r* 



,-2/3 



(86) 



where S" c is the entropy per comoving volume at t = t c and (gj) is the appropriately-averaged 
value of g* over the decay interval. In the second approximation we have introduced the 
variable u = r — t c and calculated in the limit u — > oo. Moreover we have used the relation 
Pr( mi '(t c ) = p$(t c ) and the following equation: 



S 



4 / 2tt 2 ^ 
3 V 45 



1/3 



nS/4 
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a p r 



(87) 



which is the general relation between entropy per comoving volume S and the energy density 
of radiation p r . It follows from r<j, ~ m(m/Mpi) 2 , Eqs. f!84|) and (jHSj) that the ratio of the 
entropy per comoving volume after decay to that before decay is written as 



AS = 




pi 



m 



4/3 



3/4 



(88) 



where the second approximate equality follows from Eq. ()82|) . 

If the entropy production occurs, the Universe should expand larger enough to cancel out 
the effect of the produced entropy. From this point of view, taking account of p r oc a~ 4 S 4 ^ 3 , 
we find 



(AS) 



4/3 



(89) 



where do is the present scale factor in the case with the entropy production. 

Finally, we consider the effect of the dilaton decay on the energy density of the present 
large-scale magnetic fields. We again assume that after instantaneous reheating the conduc- 
tivity immediately jumps to a large value and so the amplitude of the vector potential is 
fixed. It follows from Eqs. ©, (|54~jl — (I57j) . (|55|L and (|8~9j) that the ratio of the energy density 
of the present large-scale magnetic fields in the case the dilaton still evolves after reheating 
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and then decays into radiation, p#(L, to), to that in the case the dilaton freezes at the end 
of inflation pB{L,to) is written as 
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(90) 
(91) 

(92) 



where the last approximate equality follows from V / p<j> ps wexp(\n$-&), see Eqs. (jHJ) and 
(I39j) . In the right-hand-side of Eq. ()9()jl . the first factor is the contribution of the coupling 
with the dilaton, through which the energy density of the magnetic fields is enhanced because 
$r < 0. On the other hand, the second one is that of the produced entropy, which dilutes 
the energy density of the magnetic fields. 



C. Effect of dilaton decay 

From Eqs. ()58|) and (}9*2*|) we estimate the strength of the present large-scale magnetic 
fields in the case the dilaton still evolves after reheating and then decays into radiation. 
Figures 3 and 4 depict the magnetic field strength on IMpc scale at the present time B(t ). 
The former is for the case H m f = 10 14 GeV and m = 10 13 GeV, while the latter corresponds 
to the case H m f = 10 10 GeV and m = 10 9 GeV. In both these cases entropy is produced (i.e. 
AS > 1), and then the relation (J82j) is satisfied. The amount of the produced entropy in 
the former case is smaller than that in the latter. As described in Sec. V A, the region of 
the parameter A/c|<&rJ is constrained from the condition that the energy density of magnetic 
fields on all scales should remain much smaller than that of the dilaton before the field 
oscillation regime (tn < t < t osc ). This condition is satisfied in the region of Ak|$r| on each 
line in these figures. Moreover, Fig. 5 depicts the curves in the ifi n f — m parameter space 
on which the present magnetic fields on IMpc scale with each strength could be generated 
for the case (3 ps 5.0. In Fig. 5, the shaded area illustrates the region with m > 2H m f, where 
£r > ^osc and our analysis does not apply. In all the cases in Fig. 5, entropy is also produced. 
In this figure, we have taken the maximum of Ak|$pJ for each case. In Figs. 3—5, we have 
taken w = 0.01 and A ~ 0(1). 

The energy density of the magnetic fields is more enhanced for larger Ak|$rJ and larger 
X, while it is more diluted in the case more entropy is produced, see Eqs. (|9*0j) — (|9~2 |l . In 
fact, if (3 w 1 + Xe ps 5.0, that is, the spectrum of the magnetic fields is close to the scale- 
invariant one, and produced entropy is relatively small, e.g., in the case H- m { = 10 14 GeV 
and m = 10 13 GeV, in which AS = 4.6 x 10 6 , the generated magnetic fields on IMpc scale 
could be stronger than 10~ 10 G at the present time. Furthermore, in the case (3 ps 5.0, 
even if the entropy production factor is as large as AS ~ 10 24 , a sufficient magnitude of 
magnetic fields for the galactic dynamo scenario, B > 10~ 22 G, could be generated. More 
specifically, the amplitude of the generated magnetic field ranges from 10 _22 G to 10~ 16 G, 
namely the span favored by the dynamo scenario, if the dilaton mass lies in the range 
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9.3 x 10 6 GeV < m < 9.3 x 10 9 GeV in the case H ini = 10 14 GeV and (3 = 5 with the entropy 
production factor being AS = 5.8 x 10 15 — 5.8 x 10 24 . In the case H mi = 10 10 GeV and (3 = 5, 
the appropriate amplitude of magnetic field results for 1.7 x 10 4 GeV < m < 1.7 x 10 7 GeV 
in which AS ranges 9.2 x 10 15 - 9.2 x 10 24 . 

Moreover, as shown in Figs. 3 and 4, even if (3 < 5, that is, the spectrum of generated 
magnetic field is blue, sufficiently large seed field for the galactic dynamo scenario could also 
be generated. 

Finally, in the case X is so small {e.g., X < 10) that (3 is about unity, the spectrum is 
too blue for the large-scale magnetic field to be strong enough for the dynamo scenario to 
work. This result remain unchanged even when Ak|$r| is large and the field amplitude is 
enhanced with the relation (81) being satisfied, namely, without any entropy production. 
Thus the essential requirement is that the spectrum should not be too blue. 



VI. CASE OF POWER-LAW INFLATION 

So far we considered the evolution of electromagnetic fields in slow-roll exponential infla- 
tion models in the light of the fact that the recent data of WMAP favors models with small 
slow-roll parameters. In this section, for comparison with ^tJj we discuss it in power-law 
inflation models with the following exponential inflaton potential. 

U[4>} = Uexp(-(K(!>), (93) 

where U is a constant and ( is a dimensionless constant. In this case the spectral index of 
curvature perturbation n s is given by 

n s - 1 = -Qe v + 2-qu = ~( 2 , (94) 



1 




2^2 








^ ( 





* s ?u;- (96) 

According to the WMAP data |36j], n s > 0.93, and hence £ < 0.26. In this case the scale 
factor in the inflationary stage is given by a(t) oc t p , where p = 2/( 2 > 29. 

If power-law inflation lasts for a sufficiently long time, the dilaton will settle to the scaling 
solution |35| where U" ~ H [n f 2 with H [n f = p/t. Hence the solution of the dilaton in this 
regime is given by 

2 / Vv~\Kt\ 

$ = ^- In . 97 

Ak V P J 

From Eqs. (0) and (}9"7j) we find 

/ \ 2A/A 

f^ )= (lXM) oc a 2V(A P ). (9g) 
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Comparing Eq. (|23j) with Eq. (|9~H)l. we find 



P = £ + l. (99) 
Xp 

On the other hand, in the inflationary stage the Fourier modes of the vector potential satisfy 
the following equation: 

d 2 AAk ) r l ) 2A / -1 1\ dAAk.-n) o , . , . 

+ ^ j ' " +k 2 A I (k,7 } ) = 0, (100) 



(iry 2 A VP-1 1 ?/ dr\ 

where rj — J db/a(t) is conformal time. Comparing Eq. ()27)1 with Eq. (jlOOj) . we find 

'J^T) + 1 - (101) 

Here the solution of Eq. (jlOOj) is given by the same form as Eq. ()28j) . Since p > 1 as 
noted above, we can approximately identify Eq. (jHHJ) with Eq. (jl(Jl|) . Thus under this 
approximation we can develop the argument in the same way as in Sees. Ill— V. 

As noted in Sees. IV and V, where we considered the case of exponential inflation, 
if /3 ~ 5, that is, the spectrum of the magnetic fields is close to the scale-invariant one, 
magnetic fields on IMpc scale with the strength 10~ 10 — 10~ 9 G at the present time could 
be generated. In the case w = 0.01 and A ~ 0(1), if (3 « 5, X — A/A « 400. On the other 
hand, it follows from Eqs. ()99|) and (jl01|) that, in the case of power-law inflation, if /3 ~ 5, 
X « 2p > 58. Consequently, even in the case of power-law inflation with the maximal 
breaking of the scale invariance of primordial fluctuations within observational limit [36], 
X should be much larger than unity in order that the amplitude of the generated magnetic 
fields could take 10~ 10 — 10~ 9 G on IMpc scale at the present time. That is, the value of 
n s obtained by the WMAP data is so close to unity that the power-law inflation model 
observationally permitted is not practically very much different from slow-roll exponential 
inflation. 



VII. CONCLUSION 

In the present paper we have studied the generation of large-scale magnetic fields in 
inflationary cosmology, breaking the conformal invariance of the electromagnetic field by 
introducing a coupling with the dilaton field. First we considered the case the dilaton 
freezes at the end of inflation automatically as assumed by Ratra to see how the recent 
detailed observation of the primordial spectrum of density fluctuations in terms of CMB 
anisotropy 25], which favors slow-rollover inflation, affects Ratra's previous analysis. As 
a result we have found the resultant magnetic field could be as large as 10~ 10 — 10~ 9 G on 
IMpc scale at present for i7 inf > 10 6 GeV provided that the model parameters are so chosen 
that the spectrum of the magnetic field nearly scale-invariant or even red. 

Next we considered a more realistic case that the dilaton continues its evolution with the 
exponential potential after inflation until it is stabilized after oscillating around its potential 
minimum. It has two distinct effects on the final amplitude of the magnetic field. That 
is, the energy density of the magnetic field is enhanced as the dilaton evolves due to the 
exponential coupling, while it could produce huge amount of entropy as it decays at a late 
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time with the gravitational interaction. We have parameterized the evolution of the dilaton 
in terms of its mass, m, around the potential minimum and its amplitude at the end of 
inflation, $r, and adopted a view that it starts oscillation at t ~ m _1 , since the detailed 
shape of the dilaton potential around the minimum is not known due to the fact that the 
stabilization mechanism of the dilaton is not fully established yet, although there are a 
number of proposals. As a result we have found that the magnetic field could be as large as 
10~ 10 G even with the entropy increase factor AS ~ 10 6 provided that the scale of inflation 
is maximal and the spectrum is close to scale invariant. Furthermore the seed field for 
the dynamo mechanism could be accounted for even when AS is as large as 10 24 if model 
parameters are chosen appropriately to realize nearly scale- invariant spectrum. 

Thus the possible dilution due to huge entropy production from decaying dilaton is not 
the primary obstacle to account for the large-scale magnetic field in terms of quantum 
fluctuations generated during inflation in this dilaton electromagnetism. The more serious 
requirement is that the model parameters should be so chosen that the spectrum of generated 
magnetic field should not be too blue but close to the scale-invariant or the red one, which 
is realized only if a huge hierarchy exists between A and A, namely, X should be extremely 
larger than unity. This may make it difficult to motivate this type of model in realistic high 
energy theories. This feature is independent of whether one considers slow-roll exponential 
inflation or power-law inflation with an exponential inflaton potential as adopted by Ratra 
[l7l ] , because the latter model is hardly distinguishable from the former under the constraint 
imposed by WMAP data as far as the evolution of the dilaton is concerned. 
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FIG. 1: The curves (dotted lines and a solid line) in the H m ^ — (3 parameter space on which the 
present magnetic fields on IMpc scale with each strength could be generated (/? > 0). B(to) = 
1(T 9 G, 10~ 10 G, 1(T 16 G, 1(T 22 G (solid line), 1(T 30 G, 1(T 40 G, and 10- 50 G are shown (top down). 
The shaded area illustrates the excluded region from the observation of the anisotropy of CMB, 
Eq. (JZU). 
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FIG. 2: The curves (dotted lines and a solid line) in the H\ n f — (3 parameter space on which the 
present magnetic fields on IMpc scale with each strength could be generated (/? < 0). B(to) = 
1(T 9 G, 1(T 10 G, 1(T 16 G, 1(T 22 G (solid line), 1(T 30 G, 1(T 40 G, and l(r 50 G are shown (top down). 
The shaded area corresponds to T > 1 and illustrates the excluded region from the consistency 
condition decided by Eq. (|65[). This area includes the excluded region from the observation of the 
anisotropy of CMB. 
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FIG. 3: The magnetic field strength on IMpc scale at the present time B(to) in the case with entropy 
production. The lines are for the case H m f = 10 14 GeV and m = 10 13 GeV. ~ 1 + XXw ~ 5.0, 
(3 f» 4.5, f« 4.0, and /3 ss 3.5 are shown (top down). Here we have taken w = 0.01 and A ~ 0{1). 




FIG. 4: The magnetic field strength on IMpc scale at the present time B(to) in the case with entropy 
production. The lines are for the case _£fi n f = 10 10 GeV and m = 10 9 GeV. (3 ss 1 + Ww ps 5.0, 
(3 ps 4.5, ps 4.0, and ps 3.5 are shown (top down). Here we have taken w = 0.01 and A ~ 0(1). 



27 




4 6 8 10 12 14 



\og(H M [GeV]) 



FIG. 5: The curves (dotted lines and a solid line) in the H m f — m parameter space on which the 
present magnetic fields on IMpc scale with each strength could be generated for the case (3 ~ 5.0. 
B(t ) = 1(T 9 G, l(r 10 G, 1(T 16 G, 10~ 22 G (solid line), and 10~ 30 G are shown (top down). The 
shaded area illustrates the region with m > 2H m {, where £r > t osc and our analysis does not apply. 
Here we have taken the maximum of Ak|<£r| for each case, and w = 0.01 and A ~ 0(1). 



28 



